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1 Introduction 

The importance of conformal properties and causal structure in connection with Einstein's field 
equations (EFEs) is well known, particularly as regards asymptotic structure, see e.g. PHI E] 
and references therein. It is also well known that scale- invariant, self-similar 0], solutions act as 
important building blocks for our understanding of the asymptotic properties of non-scale-invariant 
solutions, see e.g. jS] and references therein. Indeed, the latter feature motivated the introduction 
of the so-called Hubble- normalized dynamical systems formulation of EFEs (henceforth denoted 
as UEWE), which yielded regularized field equations in the neighborhood of generic singularities 
and progress as regards our understanding, and ability to numerically handle |Hj , asymptotic 
dynamics towards such singularities. 

The purpose of this paper is to provide a framework that naturally captures all the above 
aspects. We thus combine the conformal and dynamical systems approaches to EFEs into a com- 
mon geometric instrument; conformal transformations are used to obtain scale-invariant, and thus 
dimensionless, regularized formulations of EFEs that naturally incorporate key asymptotic causal 
structures. The goal is to unravel features of the solution space and properties of solutions of 
EFEs. The relationship between conformal and asymptotic structures provides a systematic geo- 
metric foundation for finding suitable geometrically interpretable variables for dynamical systems 
analysis. This is to be contrasted with the more or less random ad hoc choices of variables that 
characterize the history of dynamical systems studies in general relativity. 

The outline of the paper is as follows: In Sec. 2 we discuss dimensionality under scale- 
transformations and use this to restrict the choice of conformal factor and frame; subsequently 
we give a set of conformal, dimensionless, field equations. As an example, we address asymptotic 
temporal structures associated with singularity formation in Sec. 3. We thus consider a timelike 
congruence and make a 1+3 split of our variables and equations. We then give some examples 
of useful conformal and temporal gauge choices: one that geometrizes and slightly modifies the 
Hubble-normalized approach used in e.g. UEWE and [S], and one that can be extended and modi- 
fied to a dimensionless autonomous first order symmetric hyperbolic system. In Sec. 4 we conclude 
with a discussion and a comparison with other work; we also outline our underlying philosophy. 

2 Scales and conformal transformations 
2.1 Scales and dimensions 

Consider a spacetime (M, g) where M is a suitably smooth 4-dimensional manifold and g is the 
physical Lorentzian metric field with signature (— , +, +, +). We use units c = 1 = 8nG so that all 
geometric properties can be dimensionally expressed in terms of a length scale. Let I be the unit 
of length, then each physical geometric field $ transforms under a scale transformation t' = SI, 
where S — const, like = S q where q defines the geometrical object's dimension, see [5] which 
we refer to for further discussion. 

General relativity is characterized by general coordinate covariance; coordinates x^ (fi = 
0, 1,2,3) are to be regarded as just labels for different spacetime events. Hence, in general, they 
do not carry any physical significance and thus it is natural to regard them as dimensionless, i.e., 
they have q = 0. This is in contrast to the spacetime interval ds 2 which describes an invariant 
physical property and naturally have a weight q = 2. Since ds 2 = dx^ dx v it follows that 
has q = 2. Hence one has to take into account the specific positioning of indices for a geometric 
object when considering its dimensional weight; for example, in the case of the energy-momentum 
tensor: T^ v , T^ v , T^ have dimensions q — —4, q — —2, q = 0, respectively. 

Although coordinates are not dimensional in general, it is natural to assign them dimensional 
weight when they express invariant physical properties that occur in special cases, which tend to 
dominate the literature. For example, the radial area coordinate when one has spherical symme- 
try; proper time or length along a given timelike or spacelike congruence, respectively; the affinc 
parameter along a geodesic congruence — these quantities all naturally carry weight q=l. Another 
important special case is weak gravity where one has a Minkowski background. This background 
is preferably expressed in Minkowski coordinates with natural weight q = 1, since they constitute 
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affine parameters of geodesies and thus express invariant features in the Minkowski spacetime. 

For computational purposes, one needs to specify a frame and use components. This introduces 
an additional element since it is quite natural to assign different dimensions to different choices of 
frames. There are three types of frames that make dimensional counting particularly easy for the 
components of geometric objects: (i) Coordinate frames, since the coordinates in general are to 
be regarded as dimensionlcss. (ii) Orthonormal frames (ONF); in this case it is natural to regard 
the constant metric coefficients as dimensionless and instead let the orthonormal one-forms carry 
dimension 1, while the dual vector fields have q — — 1. This yields that connection components 
have q = — 1 while ONF curvature components have q = —2. The ONF approach was taken as 
the starting point in UEWE, which we refer to for further discussions about dimensions in ONF 
contexts, (iii) The third choice is the one that is going to be explored in this paper — conformal 
ONF where the conformal factor carries the dimension. 

We define a conformal ONF by 

g = * 2 r ?a6 u; a u; 6 = 1' 2 g, (1) 

where r\ ab are constants (a, b — 0,1,2,3), and where g = rj ab uj a uj b is an 'unphysical' metric 
expressed in an ONF. The conformal factor ^ > 0, which is a function of the spacetime coordinates, 
is chosen so that it has dimension length, i.e. q = 1, so that the metric g (and the one-forms u>° 
and their dual frame vectors e a ) becomes dimensionless. 

Since ^ carries the dimensional scale, it follows that everything is compared with this scale 
in terms of dimensionless ratios. However, there exists no overall preferred global choice of ^ 
suitable for all possible situations; this is reminiscent to the coordinate issue — there exists no 
global coordinate choice for a general spacetime either. Precisely as in the coordinate case, one 
has to patch together a complete spacetime from regions where one has used different conformal 
factors. Instead of aiming for some global conformal choice, '5 is to be adapted to the particular 
local feature one is interested in. This can be some preferred structure associated with asymptotic 
features, or a structure associated with special initial or boundary conditions or symmetries. Useful 
candidates are obtained from suitable functions of dimensional coordinates or scalars, or quantities 
that preserve the defining structure of a special class of spacetimes. Examples are e.g. VP = i, where 
t is proper time along a timelike congruence; 'J tx where 9 is the expansion of a null or timelike 
congruence (q = — 1 for 6); = r, where r is the radial area coordinate in the case of spherical 
symmetry. We will use * = t, <J> cx associated with a timelike congruence, as examples below 
in the context of asymptotic temporal properties toward spacetime singularities. 



2.2 Conformal transformations and dimensionless field equations 

Let us begin with some notation and definitions. Consider some arbitrary metric g in some basis of 
vector fields e a , with a dual basis of one-forms {^ a }, i-e. (w a , e^) = S a b . Let us further introduce 
a connection Vfce Q = Y c a b e c , Y c a b = (lj c , V{,e a ), where V& := Ve^- We then assume that the 
connection is (i) torsion-free and (ii) metric: 

(i) V u v- V v u -[u,v] =0, (ii)Vg = 0, (2) 

where u and v are two arbitrary vectors and Vu = w a V a . The components of the metric in the 
basis {e a } are given by g a b — g(e a , e b ). Commutation functions c c a b are defined by 

[e a ,e fc ] = c c ab e c . (3) 

The curvature operator, defined by lZ(u,v) — [Vu, V<y] — Vpu .v], yields the Ricmann curvature 
tensor expressed in components as R a bcd = { <*> a , 7^-(e c , ed)e/, ) , while the Ricci curvature tensor and 
scalar, and the Einstein tensor, are defined by Rbd = Rdb = R%ad , R = R ab ab , G a b — R a b — \g a b R- 
The components of the connection and Riemann and Ricci curvature are 

Tabc = -\ [e a (g bc ) - e b (g ca ) - e c {g ab ) + c abc + c bca - c cab ] (4) 
R a bcd = 2e[ c r a | fc | d ] + 2T a f^Tl \ b \d] - F a bfC f c d (5) 
Rab = 2e[ c r C | a | fc ] + T c dc T d ab — T c a dT d bc , (6) 
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where T abc = g a d^ d bc- K follows that R a bcd = ~ Rabdc, Rabcd — ~ Rbacd, Rabcd — Rcdab, R" '[bed] — 

0, V[ e R a \b\ c d] = 0, where the two last relations are the cyclic and Bianchi identities, respectively. 
See e.g. ^U] for further discussion. 

Let us now consider a 'physical' metric g conformally related to an 'unphysicaP metric g 
according to 

£ = g ab uj a uj b = ^ 2 g ab u a uj b = * 2 g, (7) 

where uj a is a dual basis of one- forms to the basis e a : (uJ a , e b ) = 6 a b . Then the connection and 
Ricci tensor of g are related to the connection and Ricci tensor of g according to (easily derivable 
from Eqs. J3J and ©; or see e.g. 

t a bc = T a bc + 2S a (b r c) -g bc r a (8) 

Rab = Rab-2V {a r b) +2r a r b -g ab (V c r c + 2r 2 ), (9) 

in the basis e a , where 

ra--=—, (10) 

and r a = g ab r b: r 2 — g ab r a r b and V a r b — e a r b — T c ba r c . We now express the unphysical metric 
g ab in an ONF so that g ab = rj ab — diagj— 1, 1, 1, 1], and thus the physical metric is given in a 
conformal ONF: g ab = ^ 2 rj ab , see Eq. (Q. 

We take the frame variables e a M , defined by 

e a = e^d/dx" , (11) 

and the commutator functions c a bc as our basic variables, possibly supplemented by r a ; we will 
give some examples in Sec. 3. 

The governing dimensionless equations are the commutator equations, the Jacobi identities for 
e a , and EFEs, which in the conformal ONF, {e a }, are given by 

2e [a e b f = c c ab eS (12) 

e[ a C be] = C e [aC e bc] (13) 

Rab = T ab ~^ abV cd T cd , (14) 

where T ab are the dimensionless components of the energy-momentum tensor in the conformal 
ONF (recall that the energy-momentum tensor T M „ is dimensionless, i.e. g — 0, and hence it is 
possible to make the identification T ab = T ab in the conformal ONF), and where now 

Rab - Rab + U a b, U ab := -2 V (o r b) + 2r a r b - T] ab { V c r c + 2r 2 ) (15) 
r° 6c = -\v ad [Ve d c e bc + 2r le(b c e c)d ] . (16) 

Note that 2r a [ &c ] = — c abc and r a ( bc ) = — C(f, c ) a ; one can thus use T a bc as variables instead of c a bc . 

When one has a non-trivial matter source the above equations have to be supplemented with 
appropriate matter equations, however, one always have local energy-momentum conservation: 
V a f ab = 0. Recall that f» v has dimension q = -4. In the conformal ONF we have: f ab = 
~ac~bdf cd = ^-±ri ac ri bd T cd = ^- 4 T Qb ; where we have defined the dimensionless object T ab := 
rj ac r) bd T cd . This yields (see Q2]); 

e b T ab + T a db T db + T b db T ad + 2r b T ab - r a r, cd T cd = . (17) 

It is of key importance to note that the equations for the dimensional variable ^P, e a 'f = \t r a , 
decouple from the above dimensionless equations (this is seen explicitly, but also follows directly 
from dimensional reasons). The components of r a are either given functions of coordinates or 
functions of the dimensionless state space variables that depend on what type of choice of \& one 
has made, examples will be given below. This means that ^ can be obtained after one has solved 
the essential dimensionless equations, and thus ^ itself plays a 'passive' subsidiary role. It is \1/ 
that carries the scale that typically asymptotically leads to that dimensional quantities blow up. 
Since this scale now has been factored out of the problem, leaving an asymptotically regularized 
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dimensionless system, this drastically simplifies an asymptotic analysis. Once an asymptotic anal- 
ysis has been accomplished for the essential dimensionless equations, the result can be used for a 
relatively simple asymptotic analysis of the decoupled equations for thus yielding a complete 
physical result. This geometric splitting of the problem into more easily handled problems is the 
conformal ONF approach's main advantage. 

It is of interest to relate the present variables to those that one uses in the ONF approach. 
In the latter approach one uses a basis so that g = rj ab uj a uj b and an associated dual basis e a , 
(a> a , e b ) = S a b . The variables in this approach are the frame variables e Q M , defined by e a = 
Ba^d/dx 11 , and the commutator variables c a bc , defined by [e&, e c ] = c a bc e a . These variables are 
related to the present ones as follows: 

e a M = r'e/; e„'* = *e '' (18) 

c\ c = * _1 {c a bc + S a [b r c] ) ; c a bc = ^c a bc -5 a [b r c] , (19) 

where c a bc = ( ui a , [e b , e c ] ) and c a bc = ( uj a , [e b , e c ] ), i.e., the above equations are not tensor 
equations; instead they relate the components of the frame variables and the commutator functions 
of two conformally related sets of basis vector fields. The above relationships explicitly show that 
e a M and c a bc are dimensionless if VP has dimensional weight q = 1, and that everything is measured 
with respect to the scale carried by \f r . 

To make this more concrete we will consider a timelike congruence and give some examples of 
conformal and temporal gauge choices. 

3 The 1 + 3 conformally orthonormal approach 
3.1 1+3 decomposition 

We here adapt our formalism to a timelike reference congruence. We therefore choose a time 
coordinate along the congruence and align one of the basis vectors tangentially to it; this allows 
us to make a 1+3 split of the variables. 

The main application in this paper is the use of conformal regularization towards a generic 
singularity. This means that there is a close connection with UEWE. Unfortunately the notation 
in UEWE is not adapted to the conformal formalism at all, which suggests that it perhaps would be 
best to use new notation that is naturally associated with the conformal approach. Nevertheless, 
in this paper we adopt a notation that follows that of UEWE as closely as possible, since this 
emphasizes the close connection and simplifies a comparison between UEWE and the present 
work, even though this leads to some awkwardness. 

In UEWE the starting point was the ONF formalism associated with the physical metric. Since 
conformal transformations were not discussed in UEWE the frame vectors did not have hats. We 
therefore now drop the hats on the ONF vectors, i.e., {e a } — ► {e a }. This causes a problem for the 
conformal ONF vectors, however, these vectors correspond to the Hubble-normalized vectors d a 
in UEWE; we thus rename the conformal ONF vectors according to {e a } — > d a , but in contrast 
to UEWE (associated with ^ = H^ 1 where H is the Hubble variable), \& is now any function 
with q = 1. A 1+3 split (in contrast to the 3+1 split done in UEWE; see e.g. ^Hh an d references 
therein) of the ONF and conformal ONF vectors and variables yields: 

e - M~ x d t , e a = e a \M L d t +di) (20) 
do - M- X d tl 8 a =E a i (M i dt+d i ) , a = 1,2,3; i = 1,2,3, (21) 

where eo (do) is the future-directed tangent to the physical (unphysical conformal) timelike refer- 
ence congruence; M (M.) is the physical (unphysical conformal) threading lapse function, and M$ 
the dimensionless (assuming dimensionless coordinates) threading shift one-form. 

The lapse, M = A4; M. = V I' _1 M, and the dimensionless shift vector are associated with 
gauge freedom while e a l and E a l are regarded as dynamical variables, related by 



E a l = * e. 



(22) 
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The commutators are decomposed according to: 

[e ,e a ] = u a e o -[£r*a /J + 0-/ + eaM w ' r + n7 )] e ( 23 ) 
[e«,ej] = 2e Q/ 3 7 w 7 e + (2o[ a ^ 7 + n" 57 ) e 7 (24) 

[00,0a] - f/ Q d -[ft<5Q ,3 + £Q /3 + eaM^ 7 + ir)]d/3 (25) 
[d a ,df,] = 2e a ^W^d + (2A [a S^ + e aps N^)d 7 , (26) 

where the above decomposition imply the following definitions: 

(27) 
(28) 
(29) 
, (30) 

where ( ) represents trace free symmetrization. Here H = ^9 is the Hubble variable, and 9 
the expansion; a a p the shear; ii a the acceleration; cu a the rotation; fl a the Fermi rotation — all 
quantities are associated with the congruence of which e is the tangent vector field; n al3 , a a are 
spatial commutator functions, which describe the three-curvature when u a = 0; for a more detailed 
description, see e.g. [S], |13) . Analogous interpretations hold for the conformal quantities. In the 
above formulas we have adhered to the conventions used in [14j . 

Eqs. (|T§jl . P? jl -l(3"U |l yield the following relationship between the ONF and conformal ONF 
commutator function variables: 

(31) 
(32) 
(33) 

W a + R a = V(oj a + il a ) . (34) 

Let us now focus on the conformal ONF approach. From the above definitions and Eq. 11GI) it 
follows that the 1 + 3 splitted connection components of the conformal metric are given by 

Tq00 = U a , r a0 /3 = 7~LSaf3 + S Q/ 3 — tafi-fW 1 , (35) 

Tq/30 = ea/37-R 7 , r Q/ 3 7 = 2^4[ ct <5 /3 ] 7 + e 1& ^ a N s p] + ^e a psN S j . (36) 

Instead of referring to Tt, £ Q| g, U a , W a , R a , A a , N a @ as commutator function variables, one may 
refer to them as connection variables, since they describe r a b c as well as c a bc- 
The commutator equations can be written succinctly as follows: 

= {d a + U a )d Q -(Sjd -Fj)dp (37) 

= 2W a d Q -Cjd , (38) 

where 

Fj := cV - -H5j -eJ^W-* + R~<) (39) 

Cj := e Q 7/3 (d 7 - Ay) - Nj , (40) 

and this suggests that the equations can be written concisely using the above notation. 

It is natural to divide the equations into gauge equations and dynamical equations, and to fur- 
ther subdivide the latter into evolution equations and constraints (if the temporal frame derivative, 
Oq, does not appear in a dynamical equation we refer to it as a constraint equation, even though 
the spatial frame derivatives d a contain the partial time derivative dt). 

Gauge equations: 

d M a = Fj Mp + ida + U^M- 1 (41) 
= CjMp-2M- 1 W a . (42) 
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Evolution equations: 

doEj = FjEp* (43) 
d H = -H 2 -| V + IW 2 + \{d a + U a -2A a )U a -\{T m +T a a ) + \U m (44) 
= -m E Q/3 - e*\ a (2E /3)7 R s - N 0)l Ug) - 2W {a R 0) + (d {a + A {a + U (a ) U 0) 

- 3 Sa/3 + T<a/3> ~ (45) 

d Q W a = -{m8j + Fj)W + ±C a U (46) 

d A a = Fj A p -\(dp + U ){m5j + Fj) (47) 

d Q N af} = -(3n5^ a + 2F^ a )N^ + (d 1 + U 1 )e^ a F s ' 3) . (48) 
Constraint equations: 

= CjEp 1 (49) 

= m 2 - iS a/3 E a " + W 2 -2W a R a + ± 3 K-T QQ + ±(U o + U a a ) (50) 
= -2d a H + {d -3A )Hj + ej~< (£/ N s ^ + 2U W 7 ) + Cj W -T 0a + U 0a (51) 

= (dp - 2Ap) Nj + eJWpA^ - 2Fj Wp (52) 

= (d a -U a -2A a )W a , (53) 

where 

3 S a p = d {ot A 0) -(d 1 -2A 1 )N s{a e 0) "< 5 + B {a0) (54) 

3 TZ = 4d a A a - 6A 2 - \B a a (55) 

B a p = 2NJ N 7 p - N a p (56) 

Uoo = -3(d +H)r a + ^(dp + 2rp) + 3LP - 2^] r 7 (57) 

U 0a = -2{d a -r a )r Q +2(mJ + + e a ^W^)rp (58) 

f^(Q/3) = 2 [£ Q(3 r - + A( a )rjj) - e> S ( a N 5 p)^ + r^rp)] (59) 

U 00 + U a a = 6(m + r o )r o -2[S " / (2Bp+r )-4A^]r^, (60) 

where we have used the notation v a v a = v 2 . If M a = = W a , then 3 TZ and 3 iS Q /3 are the curvature 
scalar and trace free part of the Ricci tensor, respectively, of the conformal 3-metric. 

A conformal 1+3 split of the equations for T ab yields (to avoid clashes with the notation in 
UEWE we will refrain from explicitly splitting T ab in terms of its irreducible parts, but see the 
discussion below): 

{do + m)T OQ +HT' 3 fj -(d p + 2U l3 ~2Ap)T ' 3 + E^ 7 T^ + C = (61) 
(0 O + 4W) T 0a - U a Too ~ A a T% + Zj T ofj - (dp + Up- 3Ap) Tj (62) 

+ e a ^[N 5 T Sl -(W 1 -R y )T o0 ]+C a = 0, (63) 

where 

Co := (T OQ + T )r o -2T o rp (64) 
C a := 2T 0a r - (T 00 - T 13 p) r a - 2Tj r . (65) 
As done in e.g. UEWE, the energy-momentum tensor can be 1+3 splitted according to 

f ab = pu a u b + 2u a q b) +ph ab + 7r ab , (66) 

where u a h ab — 0, u a q a — 0, 7r a a =0. In and UEWE the following normalization was intro- 
duced: (f2, P. Q a , II Q ^) = (p,p,q a ,ir al3 )/(3H 2 ); the reason for this is that this yields the standard 
definition of the important cosmological dimensionless density parameter f2. However, from a con- 
formal geometric perspective it follows that if one wants to use as conformal factor, then 
the natural normalization factor is H~ 2 . This then suggests the following new definitions in the 
Hubble-normalization case: (D,P,Q a ,H a ^) = (p,p,q a ,n a ^)/H 2 , associated with the irreducible 
1+3 decomposition of T ab , and hence f2 = D/3; alternatively one may use (-/^H)^ 1 as conformal 
factor, but that changes the conventions with respect to |S] and UEWE as regards the connection 
variables. 
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3.2 Examples of conformal and gauge choices 

Associated with a choice of a conformal normalization factor 'J, there exists a natural temporal 
gauge choice — the conformal 'proper time gauge', A4 — 1, which in the 3+1 case Mi = = W a 
reduces to the conformal Gauss gauge (see [3] for a discussion about the use of conformal Gauss 
coordinates to cover large spacetime domains), however, it is of course not necessary to choose 
this gauge. Instead of a general discussion about conformal, frame, and coordinate freedom, we 
will focus on a few examples. Since the emphasis in this paper is on the conformal approach, we 
will divide our discussion in terms of conformal choices; we consider two such choices — conformal 
Hubble-normalization and conformal proper time normalization. 

3.2.1 Conformal Hubble-normalization 

The first example is given by 

* = H~ x , (67) 

where H = \6 is the physical Hubble variable and 9 is the physical expansion, defined by 9 = V a u a , 
where u — eo (in the notation of Sec. 2). In this case we define the physical deceleration parameter 
q (not to be confused with the scale weight q) and logarithmic spatial frame derivative r a by 

d H = -(l + q)H, d a H = -r a H, (68) 

i.e., in terms of r a we have: tq = 1 + q = —daH/H and r a = —d a H/H, which combined with Eqs. 
(jnZJl-ljlini) determine U ab . 

Eqs. (|31(l - (|34|l and the above definitions yield: 

H = -q , U a = U" -r ai A a = A% + r a , (69) 

where := u a /H, A^ := a a /H, while the other variables are just the usual Hubble-normalized 
variables used in e.g. UEWE. 

One can choose to let q and r a be determined by the Raychadhuri equation (|44(l (the time 
derivative of q drops out when TL = — q, ro = l + q are inserted in 1)44(1 1 and the Codacci constraint 
(|5T)l (the spatial derivatives of q drop out), respectively. However, since — q is just one of the 
connection variables in the conformal formulation, it is quite natural to extend the normalized 
state space to include q and r a (which is also connected to the present formalism through its link 
to the gauge quantity U a ) as independent variables, something which has been found to be quite 
useful, see e.g. |Hj and |T5] . 

Setting Ad = 1 yields the separable volume gauge, see UEWE, and if one in addition sets 
M a = = W a one obtains the inverse mean curvature gauge, which in the present context can 
be interpreted as the conformal Gauss gauge associated with W = H , something which is also 
reflected in that the congruence is conformally geodesic: U a = 0. This further emphasizes the 
geometric nature of the present approach, and the preference of using U a instead of . 

Note that with the current choice of conformal factor, and a conformal Gauss coordinate choice, 
the present formulation reduces to the 'standard' Hubble-normalized formulation in the spatially 
homogeneous (SH) case, since U a — = r Q , i.e., the present formulation yields a natural geometric 
generalization of the Hubble-normalized SH case, discussed extensively in e.g. [5]- For the same 
reason the present approach reduces to that used in e.g. UEWE for the SH part of the so- 
called silent boundary, where the attractor for generic singularities resides (this is also the case 
for the subset associated with isotropic singularities |16p. Hence the description of the attractor 
for a generic singularity in the present geometric formulation is identical to that in UEWE — the 
asymptotic regularization properties are generically identical. 

3.2.2 Conformal proper time normalization 



The second example uses the physical proper time i along a timelike reference congruence as the 
conformal factor: 

* = i\ (70) 
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The time variable is subsequently reparametrized so that one obtains a dimensionless time variable, 
t, according to 

t = ln(t/t ), i=t e\ (71) 

where io is some reference time; it follows that the new time variable is just the conformal proper 
time, since M — 1, and hence So = d/d t and r a = (1, 0, 0, 0). This leads to 

C/oo = -m , U 0a = 2U a , U {aP) = 2S Q/3 , C/oo + U a a = &(2H + 1) , (72) 

and thus one obtains a first order autonomous system of equations. This can be seen explicitly, 
but again this also follows from dimensional reasons: t is the only varying quantity that carries 
dimension and t can hence not appear in the dimensionless equations, and therefore the same 
holds for t; neither does the normalization affect the essential first order structure of the usual 
dimensional ONF approach. 

In ^7], [T^l, and it was shown that by extending the ONF approach to also include the 
curvature tensor and the Bianchi identities one can obtain a first order symmetric hyperbolic 
system, if one uses proper time along a timelike congruence (this was shown for a perfect fluid 
with a barotropic equation of state by using proper time along the fluid congruence). It is of 
course also possible to extend the current conformal ONF approach similarly, as will be discussed 
in the next section. Since if? = i = t^e 1 does not modify the principle parts of the equations of 
a curvature extended formulation, we conclude that it is possible to extend the present 'minimal' 
formalism and obtain a dimensionless autonomous first order symmetric hyperbolic system for the 
i-normalized equations; incidentally, this system is of course well-posed. 

Let us now for simplicity specialize the temporal reference congruence to be non-rotating, 
M a = = W a , so that we obtain proper time normalized equations and a conformal Gauss 
coordinate system, for which U a = and hence Uo a = 0. Moreover, let us consider a generic initial 
spacelike singularity and let us specialize the time coordinate so that it becomes a simultaneous 
bang function, i.e., t = at big bang, and hence t — > — oo towards the singularity. We thus 
take the synchronous coordinates used by Belinskii, Khalatnikov, and Lifshitz (BKL) [20] as the 
starting point (see also |2J for a discussion about the existence of such coordinates) and obtain a 
dimensionless formulation that brings us particularly close to the work of BKL, which therefore can 
be interpreted directly in terms of the dimensionless state space picture the present formulation 
gives rise to. 

It may seem that the latest approach is superior to the Hubble-normalized one, however, both 
have advantages and disadvantages. The advantage of the Hubble-normalized approach is that one 
essentially uses the expansion which appears prominently in the singularity theorems, and that 
one decouples the dimensional variable H. This implies that H carries the dimensional constant 
of integration, which we denote as the scale parameter even though it is a function in general, 
when one has a scale-invariant source. The advantages of the conformal proper time normalization 
approach is that it yields a first order system which may be extended to a first order symmetric 
hyperbolic system, and that one obtains a formulation closely related to that of BKL, which 
facilitates comparisons. A disadvantage is that one does not decouple a variable that carries the 
scale parameter. The difference of the two approaches as regards the last aspect can be illustrated 
by the Kasner subset. 

Let us for simplicity only consider the vacuum case (as discussed in UEWE, if one has a source 
one may have additional test fields such as the 3- velocity of a fluid). In the Hubble-normalized 
approach the Kasner subset is defined by setting all variables to zero except the shear which satisfies 
1 — gS Q( 3 E"' 3 = (and q = 2 if we consider the r a extension), which yields the so-called Kasner 
sphere, see UEWE. On the other hand, when we use the conformal proper time normalization, 
with a simultaneous bang function, then Ti and T, a p ^ 0. Setting r a = (1,0,0,0) and all other 
variables to zero, apart from Ti and S Q| g, leads to that Eq. (j50|l yields: (1 + Ti) 2 — g£ a /3 £ = 0, 
i.e., we obtain a cone with Ti = —1 as apex. If we consider expanding models, then Ti > 0; for 
Kasner Ti — Tio — const > 0, and hence we obtain a Kasner sphere for each value of Hq. This 
illustrates that in contrast to the Hubble-normalized formulation, we obtain a scale parameter as 
a constant of integration in the proper time normalized formulation, something that somewhat 
complicates the description of the structure of the attractor for generic singularities. 
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Implicit in the above discussion is also the need for choosing a 'dominant' quantity as conformal 
factor in order to obtain asymptotically regular and well behaved equations, e.g., if we had used 
the inverse of a component of n a ' s as the conformal factor many state space variables would have 
blown up towards an approach to Kasner (and towards a generic singularity). Thus, e.g., for a 
generic singularity one needs to use a conformal factor that goes to zero at least as fast as H^ 1 
in order for the state space variables to remain finite, however, it is preferable if \& H remains 
finite towards the singularity, i.e., it is preferable to have a 'marginally dominant' conformal factor 
that leads to finite state space variables, without all of them going to zero, and well behaved field 
equations. 

4 Discussion 

In this paper we have used conformal transformations in order to obtain dimensionless regularized 
field equations that allow one to extract asymptotic features and properties about the solution 
space of general relativity. The conformal factor is to be chosen so that it captures a characteristic 
scale associated with asymptotic structure so that all the state space variables form dimensionless 
ratios with respect to this scale. In this paper we have used a 'minimal' approach, however, we 
could have extended our formalism to also include the curvature, in particular the Weyl curvature, 
and the Bianchi identities, as done in e.g. ^3 and ^Hj- Since the Weyl curvature is conformally 
invariant, this implies that one should use the Weyl curvature in a conformal ONF, i.e., in contrast 
to Friedrich's conformal approach, see |TZ] and references therein, the Weyl tensor is not to be 
scaled with the conformal factor in our approach — it suffices to express it in a conformal ONF. We 
here note that even though the conformal factor enters our equations implicitly in the combination 
e a ^ = r a (Sec. 2 notation), the components of r a stay finite when $ — * 0, if \& is chosen as 
an appropriate marginally dominant scalar that carries dimension q — 1. This leads to a coupled 
system of regular dimensionless field equations, since the equations (differential or algebraically 
trivial, depending on the choice of VP) for the dimensional '5 decouple; furthermore, it is the reduced 
dimensionless system that carries the essential dynamics, since "J can be obtained afterwards once 
the equations of the reduced system have been solved. 

Note that the components of the Weyl tensor in an ONF (which have dimension q = —2) and 
the components of the Weyl tensor in a conformal ONF (dimension q = when the conformal 
factor have dimension q = 1) only differ by the square of the conformal factor. In and [2] 
the Weyl tensor was normalized with (-\fSH)~ 2 as the scale factor, since the same factor was used 
to normalize T ab , as discussed previously. However, we now see that from a conformal point of 
view the natural normalization factor is just the square of the conformal factor, which in the 
Hubble-normalization case is H~ 2 . 

Our choice of conformal factor is also quite different than that used in conformal approaches 
to isotropic singularities, see e.g. [22], [S]. There the motivation for the conformal factor is a 
purely mathematical one; choose a conformal factor so that regular expressions for the covariant 
coordinate components of the 3-metric at the singularity are obtained. This typically leads to 
a conformal factor that does not have any particular dimension, indeed, the dimensional weight 
is different for different matter sources (not surprisingly, increasingly complicated dimensional 
conformal properties lead to increasingly messy subsequent mathematical analysis — this is why, 
e.g., the dimcnsionally simple case of a perfect fluid with radiation as equation of state is relatively 
easy to treat). The present approach uses a strategy that is almost the opposite. The conformal 
factor is always chosen to carry the dimensional weight and for isotropic singularities, see as 
well as for typical timelines for generic singularities, all the components of the covariant 3-metric 
blow up, and this is a very good thing! Instead the focus is on the components of the spatial frame 
vectors which determine the contravariant components of the 3-metric; these components all go to 
zero, and this directly reflects the asymptotic causal properties towards the singularity — asymptotic 
silence, see UEWE and 7 . The present conformal approach emphasizes the geometrical content 
of the discussion about asymptotic silence in UEWE and even further due to the connection 
between causal and conformal properties. In the present approach the focus is on the conformal 
state space which is extended to include the so-called silent boundary where all components of the 
contravariant 3-metric are zero. This extension then allows one to use the state space picture to 
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perturb the structure on the silent boundary into the physical state space, and thus derive physical 
results about asymptotic spacetime properties. 

In this paper we have used temporal asymptotic aspects associated with singularities as an 
example and made contact with other work to illustrate the usefulness of our approach. However, 
it is our belief that the current formalism may be useful for all types of asymptotics in general 
relativity: temporal, null, and spacelike; for non-isolated and isolated systems. Indeed, we already 
know that it is useful for future temporal asymptotes in SH contexts, since it contains the Hubble- 
normalizcd formalism which already have proven to be useful in this regard. However, it should be 
pointed out that the Hubblc-normalization did not lead directly to regularized equations towards 
the future for the general SH models, some additional manipulations were needed in order to 
obtain asymptotic results, but the Hubble-normalization provided the first key step |24| . A similar 
situation is expected for null infinity. The present formalism is expected to yield direct results for 
special cases, but not the most general ones where additional manipulations will be necessary. The 
current formalism could have been used as the starting point in the work |25| , using the radial area 
coordinate r as conformal factor, which yielded asymptotic results as regards spacelike asymptotes, 
for non-isolated and isolated systems, in the context of static spherically symmetric spacetimes (in 
[23] there existed two relevant scales, however, several scales can be handled by first using the 
conformal transformation to take care of an overall scale, and then making additional variable 
transformations that form further ratios, which one by one quotient out the other scales). 

However, the main reason for believing that our proposed approach may be a useful ingredient 
in future studies of asymptotics perhaps comes from the simplicity and naturalness of the main 
ideas — summarized as follows: 

(i) Consider a marginally dominant dimensional scale that captures some key asymptotic fea- 
tures. 

(ii) Use conformal transformations to geometrically quotient out and decouple this scale so that 
all remaining quantities represent dimensionless ratios with respect to that scale. 

(iii) Use the obtained reduced dimensionless regularized field equations on an extended state space 
(i.e., include asymptotic limits if they occur on the boundary of the original dimensionless 
state space) to derive and describe asymptotic properties. 

Thus dynamical systems approaches in general relativity, based on regularized dimensionless field 
equations, have found their place in a familiar conformal geometric setting. 
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